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Introduction 

The integrable structure of the Calogero-type systems is of special interest concerning 
such diverse problems as hierarchies of partial differential equations 0, Seiberg-Witten 
solutions of N=2 SUSY [§, quasi-Hopf deformations of quantum groups and many 
others. There are many different approaches to describe the Hamiltonian structure of 
such systems. One of the aims of this paper is to combine the Hitchin description and the 
language of the universal algebraic-geometric symplectic form. The method of the special 
inverse spectral problem was successfully applied in the matrix KP case in 0. This 
method deals with the auxiliary linear problem coming from the Lax representation of 
the equation. The Baker- Akhiezer function (FBA) turns to be a function on the spectral 
algebraic curve, and its specific analytic properties are very restrictive and provide the 
reconstruction of the potential and the poles of FBA in terms of theta-functions on the 
Jacobian of the spectral curve. The same procedure was applied to the spin Ruijsenaars- 
Schneider system in ||^ . But for a wide class of integrable systems of the Hitchin type the 
explicit solutions remain still unknown. In the first section we revisit the construction of 
the Hitchin system on the degenerate curve and especially the Gaudin elliptic two-point 
N-particle system with spin. In the second section we investigate the finite gap potentials 
of the matrix Davey-Stewartson equation. The system of non-linear equations of its poles 
dynamics is deduced and the algebraic-geometrical correspondence between the phase 
space of this system and the space of spectral data is constructed. In the third part the 
FBA is reconstructed in terms of ^-functions, the explicit expressions for the potential 
and for the solution of the poles dynamics are given. The universal symplectic structure is 
constructed in the fourth section. We also demonstrate that the system of poles dynamics 
is hamiltonian and identify it with the Gaudin elliptic system with spin. In the last section 
we discuss the rational and special modular degenerations of the system. 



1 Gaudin system 

The Hitchin system was introduced in as an integrable system on the cotangent bundle 
of the moduli space of stable holomorphic bundles on the algebraic curve S. The phase 
space was obtained by the Hamiltonian reduction from the space of pairs (i^, $, where 
dj^ is the operator defining the holomorphic structure on the bundle V and $ is the 
endomorphism of this bundle, more precisely $ G ^"'^(2, EndiV)). The reduction by the 
action of the loop group GLj^{z) was done, i.e. the group of GLjv- valued functions on S. 
The invariant symplectic structure on the base space writes 



The zero level of the moment map is described by the condition = which means 
that $ is holomorphic with respect to the induced holomorphic structure on the bundle 
Endiy). It turns out that the system of quantities Tr^'', treated as vector functions on 
the phase space, Poisson-commute and their number is exactly half the dimension of the 




(1.1) 
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phase space. This system was stated in a formal way and its explicit formulation was 
obtained for special cases in P, pUj . 



The generalized construction of the Hitchin systems on degenerate curves [||] provides 
a more explicit description. One considers the moduli space of semi-stable bundles V on 
an algebraic curve with singularities. The Higgs field in this case is a meromorphic section 
of the bundle End{V). The holomorphic bundle on the singular curve is the holomorphic 
bundle on its normalization and the set of additional data characterizing the bundle over 
singularities. 

Let us consider the elliptic curve with two simple fixed points zi,Z2. The Higgs field is 
a meromorphic section with two poles on the curve. The moduli space Ai of semi-stable 
bundles on the elliptic curve with such type of singularity decomposes into the moduli 
space of semi-stable bundle V on the normalized curve and the finite dimensional space 
of additional data at the fixed points: Qi G End{V\zi)- The cotangent bundle T*M. 
has additional coordinates pi e T* End{y\zi) where i = 1,2. The modified symplectic 
structure is 

y 5$ A + Tr{6{g^'p^) A 6g, + 6{g^'p2) A 692) (1.2) 

and the zero level of the modified moment map is 

9$ + [A, $] + pi6\z,) + p2S\z2) = 0. (1.3) 

We restrict ourselves to the subspace of semi-stable holomorphic bundles of rank which 
decompose into the sum of line bundles. In the basis associated to the decomposition we 
can make the matrix A diagonal with elements Xi, . . . , x„. The equation (|1.3| ) becomes 

d<l>ij + {xi - Xj)<l>ij + {pi)ij6\zi) + ip2)ijS\z2) = 0. (1.4) 

It could be solved as in |jl| in the following way 

exp{{xi - Xj)f5f ) a{xi - xj + z - zi) a{xi - xj + z - Z2) . 

= {{Pl)ij^ Vt V + (P2)ii^ Vt ^^-^^ 

zTTi a[Xi — Xj)a{z — zi) a[Xi — Xj)a[z — Z2) 

^ii = Ui + {pi)iiC{z - zi) + {p2)iiC{z - Z2); (1.6) 

where r is the modular parameter of the elliptic curve, a is the WeierstraB cr-function and 
Ui are additional free parameters. Now res2=z-^47r^Tr$^ is equal to 

Hi = \2 UJi{pi)ii + y'(pi)n(p2)nC(2;i - Z2) + y](pl)ii(p2)ji^T^^^ Xj + Zi Z2)^ ^^^^^ 

and the residue at the point z = Z2 gives a similar expression for the Hamiltonian H2 
which differs from Hi by the interchanging of zi and Z2- These functions are called the 
Hamiltonians of the Gaudin elliptic system with spin. We are interested in the complete 
integrability of this system. On the moduli space of holomorphic bundles on the regular 
curve in |^ the number of coordinates on the set of quantities Tr^^, i. e. the number 
of Beltrami differentials, was proved to be exactly half the dimension of TA4. Our case 
differs from the classical one and we calculate the dimension directly. The coefficients 
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of the principal part of Tr^'' at the vicinity of the fixed points are the integrals. The 
mentioned quantities are elliptic functions, so the sum of their residues is equal to zero. 
The number of integrals is equal to N{N + 1). The number of coordinates on the phase 
space Xi,uJi{pi)ij, {p2)ij, {gi)ij, {92)ij is equal to 2iV + 4A^^. The symplectic structure ( |1.2| ) 
is invariant under the action of the group GL^^ x GL^^ in the following way: the element 
(91,92) acts as 

9i ^ 9i9i, 92 ^ 9292- 

The moment map of this action is ^ = {{9i)~^Pi9i, {92)~^P292) and we can fix it to have 
a diagonal form and this kills 2N{N — 1) degrees of freedom. The stabilizer of this choice 
of moment map consists of the diagonal subgroup Diag x Diag E G x G. It means that 
the choice of the representative of the stabilizer's orbit in the moment map subspace kills 
2N additional degrees of freedom. The dimension of the reduced phase space is equal to 
2N + 4A^2 _ 2N{N - 1) - 2iV = 2A^(A^ + 1) so that the number of integrals is exactly half 
the number of variables. 



2 Davey- Stewart son equation 

As was shown in |TT|, the Davey-Stewartson equation 



iUt + Ua:x - Uyy + U(j) = 4>xx + 4>yy = '2{\u\l^ - \u\l ) (2.1) 



has a Lax representation L = [L, M] with 

^=(0 0) ' = ^ + 'y- 

Let us consider the /-dimensional generalization of the auxiliary linear problem related to 
this operator 

where Ai,A2 are / x /-matrices, di = E,d2 = —E, and E is the identity / x /-matrix, 
and \1/ is a 2/-dimensional vector function. Looking for the quasi-periodic solutions of the 
matrix Davey-Stewartson equation we put the potential in the special form: 

-"■"I?"'*' .(',Mx-x,)- ''^--l^oA _,^.). (2.3) 

where a;[ and bl. are /-vectors, and a is a WeierstraB function. As usual we look for the 
quasi-periodic solutions of the linear problem. The convenient basis of such functions 
with simple poles on x is provided by 

^i(z,x) = ""1^, ^\ exp{-({z)x), where Z2 - zi = rj. 

(7[x)a[z - Zi) 

We call the first / components of the function and \l/2 the last / components. The 
blocks of the potential have a monodromy such that the ratio of the monodromies of 



3 



and \l/2 must be equal to the monodromy of Ai. We look for the eigenfunction in the 
form: 

n ^ 

*i = Sl<^i{z, X - Xj) exp{ky+ + -y_), (2.4) 

where y± = ^{t ± x). We also look for the solution of the dual equation 

^^J+DA = (2.5) 

in the form 

"^t = Yl ^tMz, Xj - x) exp{~ky+ - -yJ). (2.6) 



Theorem 1 The equation ( ^.^ has a solution of the form l \2.4\ ) and the equation 



has a solution of the form l \2.(^ if and only if the poles dynamics is described by the 
system: 

a{ + A{a{ + Y: a[{b[\-)^^l^f-^^^ = 0; (2.7) 
^ a{r])a{xj - xi) 



,+ i.cr{xj ~xi + r]) 



Xj 



cr{r])a{xj - xi) 

^-^^^T^^^^ t^;^:;!:!, -^^ (2.9) 

air])aix, - xi) ^' a{r])a{x, - xi) ) ' 

K - K = E(^2^al)(^'1^4)^||^=^^(C(^ - xi) - C{x, -x, + r^) + 2C(r/)) + 

|J(6i\^)(6r«l)^||^=^^(C(x, - - C(x, -x,-v)- 2C(r/)). (2.12) 

Proof As was done in [P, ^ |l^], introducing the expression (p.4|) into equation (|2.2| ) and 
eliminating the essential singularity, we calculate the coefficients of the expansion in x at 
the vicinity of the poles x = Xj. The lowest term is of order —2 

Siixj - 1) + a{b{^ si = 0; 

Si{xj + l) + aibCsi = 0. 
The similar condition from the dual equation gives 

st{x,-l) + staibC = 0; 
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si {xj + 1) + si a{b{^ = 0. 
We can fix the gauge of variables a;^, 6;^ such that 



1-Xj = (6{ ai), -1 



X 



and so 

Si = <'oi. sC = d% k = 3-k, (2.13) 
where are scalar functions. The terms of order —1 at x = Xj implying ( p. 131 ) can be 



represented as (|2.7|) , (|2.^ ) where A;^ is defined from the following equations on c* 

{dt + kl + Li)C = 0; (2.14) 

{dt + ^I + L2)C = 0, (2.15) 

where 

(Li),,. = 5,,((6i+4)(C(^ - z,) - C{z) - Civ)) - Al) + (1 - 6,,){b\^ai)^,{x, - x,, z); 

iL2\^ = 6,M''<)iCiz - z,) - C(z) + av)) - + (1 - 6,,){U,^a{)^,{x, - x„ z). 

Concerning the dual equation, we obtain similar equations for the variables and the 
equations (p.9|) , (|2.1CI|) . Thus for the Lax equation we have 

9t(Li-L2) = [Li,L2], 

which is equivalent to the equations (|2.11|) ,( pn^ ). ■ 

The standard trick in the method of algebraic-geometrical direct and inverse problem 
is that the spectral curve F, defined by the characteristic equation det{k — ^ + Li — L2) = 0, 
has specific properties and that \E^-function defines a holomorphic line bundle on the curve, 
and its essential singularity depending on time is the deformation of the fixed line bundle. 
The analysis is rather traditional and we venture to outline main propositions on the 
analytic properties of eigenfunction and properties of the spectral curve. We use the 
notation L{z) = Li{z) — L2{z) and k = k — First, noting that the matrix reSz=ZiL 



z 



differs from the unity matrix by a matrix of rank / we find the decomposition at the points 

^1, Z2 



N 



RCk, z) = + ^- + hl'\z)) (2.16) 

i=l 

where h]''^{z) are regular and ul''^ = 1 if i > I. Using that F is an A^-sheeted covering of 
the base curve z and the definition of the branch points of the curve, we count the number 
of branch points; it is equal to 4NI — 2l{l + 1), and by Riemann-Hurwitz theorem the 
genus of F is = 2NI — / (Z + 1) + 1. In virtue of the Lax representation the parameters of 
the spectral curve are integrals of motion. If, in addition, we fix c{t)\t=Q = const it can be 
deduced that the poles of vector c do not depend on time. In general position the number 
of such poles is half the number of branch points, and it is equal to 2NI — + = g — 1. 
The function has essential singularities at the points over z = zi, z = Z2 and additional 
21 poles. Summarizing the analytic properties of the eigenfunction we obtain 
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Theorem 2 The function \1/ is meromorphic on the curve T of genus g outside 21 points 
Pi, and has g + 21 — 1 poles 71, 7c,+2/-i which do not depend on t,x. At the points 
Pi,i = 1, over z = z\ it has the following decomposition 

00 

^a{x,t,P) = (x?^ + ^xr(^,^)^')exp(A,(x + i))^i(0,0,P) 

s=l 

and when i — I + 1, ...,2l at the points over z — zi we have 

00 

*„(x,t,P) = (x^^ + J]xf(^,^)-2')exp(A,(x-t))*i(0,0,P) 

s=l 

where \i{z) = (1 — Vi)z~'^ — hi{z), hi are regular functions in the vicinity of z — 0, Vi 
are eigenvalues of the matrix L at the points over z — Z\^ z — Z2 respectively and Xo' ^'^^ 
constants. 

The similar proposition on the analytic properties of the solution of dual equation is true 

Theorem 3 The solution of the dual equation \E'"'' is a meromorphic function on the same 
curve r of genus g outside the points Pi, i = 1, ...,2l It has g + 21 — 1 poles 7^, 7^2«-i 
which do not depend on t, x. At the points Pi,i = 1, / over z = zi it has the following 
decomposition 

00 

t, P) = (xo+'"^ + ^^'"(^' exp(-A,(a; + t))vl>+(0, 0, P), 

s=l 

and at the points Pi over z — Z2, i — I 21 one has 

00 

*+(x,i,P) = (xo^'"^ + ^x:'"^(^,^)^^)exp(-A,(x-i))*+(0,0,P), 
where Xo'°'^ '^f^ constants. 

3 Reconstruction formulas 

In this section we demonstrate that the analytic properties that we found in the previous 
section are sufficiently restrictive and allow us to reconstruct the eigenfunction in terms of 
the characteristics of the curve, such as abehan differentials, ^-functions and Abel trans- 
form. The standard proposition about the uniqueness of the function with the following 
properties 

1. '^ai'X-, t: P) has g + 21 — 1 poles 71, 7g+2i-i on the curve F of genus g 

2. it has essential singularities at the points Pj,j = 1, .., / of the form 

^^{x, t, P) = e-"(-+*)(5«,- + P4.(x, t)w^) 
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and at the points Pj,j = / + 1, 21 of the form 

s 

where Wj{P) is the local parameter at Pj 

arises from the Riemann-Roch theorem. And its existence follows from the construction: 
let dQ} , dVt^ be unique meromorphic differentials, holomorphic outside the points Pj such 
that: 

dVt^ = d(w~^ + 0{wj)) at the points Pj 

d{wj^ + 0{wj)) at the points Pj with j = 1, / 
d[—wj^ + 0{wj)) at the points Pj with j = / + 1, 21 

normalized by the condition J^ dQ^ = 0. We set Uj''^ = ^ J^^ dQ^; then, up to a 
normalization constant 



dQ^ 



where v4(P) is the Abel transform, K is the vector of Riemann constants, 

3-1 9-1 
g+2Z-l 2« 

Zo = -K- 5^ A(7,) + 5^A(P,), Z^ = Z,-A{P^), 
i=i j=i 

Q,^'^^ = jp^^ dVt^ . The dual eigenfunction is constructed analogously 
where 



(3.1) 



g(A(p) - [/(^x - + zt)e{z^) n,^. g(A(p) + 



9-1 

5+ = -ir - A(P,) - 5^ A(7+), P+ = -K- A(7;_i+,), 

K,= 5^ (A(7.) + A(7+))-25^A(P,), Zo+ = Zo-2ir-iro, Z+ = Zo+-A(P,). 

j=i j=i 

These functions satisfy the auxiliary linear problems ( |2.2| ), ( |2.5| ). The finite gap potential 
of the Dirac operator, i.e. the solution of the matrix Davey-Stewartson equation is given by 
A = aF^a—F^, where a is the diagonal matrix with diagonal elements —1, . . . , — 1 

I I 

and Pi is the matrix coefficient of the decomposition near the singularities of the FBA. 
Therefore the solution of the system of poles dynamics ( [^.7| ),...,( |2.12| ) is given by the 
formulas: 
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1. the poles of the eigenfunction can be found as the solutions of the equation 

e{u'^^^Xi{t) + u^^h + Zq) = 0, 

2. the spin variables are given by: 

^' 9iZ+)9{S+) 

3. the normalizing factor is: 



21 



a=l 



2 ^ eiz^)eiz+)eis^)eisi) 

In the formulas above we used the notation: 

4 Universal symplectic structure 

The systems obtained as poles dynamics in the finite gap theory were often introduced 
from physical motivations. Such an explanation leaves the structural questions concern- 
ing dynamical systems apart. By structure we mean the algebraic symmetries in the 
broad sense. The importance of the hamiltonian structure cannot be overestimated in the 
fundamental problem of quantization, likewise the "action-angle" variables have a crucial 
meaning in the Seiberg-Witten theory. Heuristically most of the systems involved are 



hamiltonian and moreover integrable in the sense of Liouville. In [14, 15] was proposed 
the universal hamiltonian description of such systems in terms of the flows on the space of 
quasi-periodic operators generated by the Lax equation. First, this approach was applied 
to the description of the hamiltonian structure of the nonlinear equations of the soliton 
theory, but in [|l^] the efficiency of this method when applied to the finite-dimensional case 



was confirmed. The essential point is the algebraic-geometric correspondence which maps 
the coordinates of the phase space of the system of poles dynamics to the space of spectral 
data, namely to the set of modular coordinates, the divisor of poles of the Baker- Akhiezer 
function and the type of singularities, which are given by a pair of Abelian differentials 
on the curve, in our case by the differentials dk, dz. 

The coordinates on the phase space of the system (2.7),...,( 2.12| ) given by 



Xj, Aj, a{, Og, 6i, &^ are invariant under the following change of variables 
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where Wk, k = 1,2 are constant matrices. We also have the symmetry of the system 

where fXi are constant scalars. We also recall the condition {b{^a{) — {b2^a{) = 2. The 
dimension of the reduced space is 4M + 2A^ - 21(1 - 1) - N - N = ANl - 2l{l - 1). We 
fix Wi , W2 such that the matrices 

46*2^1^1, W^^ 4b\^W2 (4.1) 

i i 

are diagonal and the matrices Wk leave the vector (!,...,!) invariant. Further we use the 
notation A^, BJ^ for the new variables. The normalization (|4.1| ) takes the form: 

i 
i 

For the eigenvector of the operator L the following representation is true: 

Ci = {z- zi){4 + 0{z - zi))e^\>{-XiC{z)), if j < /, 

Ci = {z- Z2){4 + 0{z - Z2))ey.Y>{-XiC{z)), if 21 > j > I, 
c, = {ci + 0{z))exp{-x,C{z)), if J > 21. 

The vectors are proportional to the vectors . We chose for normalizing the following 
condition 

N N 

J2 Aimci = Sj^, J < 21, Yl ^irnci = 0,J> 21. (4.2) 

j=l 4 = 1 

It implies the customary normalizing of FBA at the points Pj 

*g = {5jq + 0{z))exp{kjx). 

Now we define the "action" variables as the poles of the vector C which solves the equation 
{k + L)C = and satisfies the normalizing condition 

IN 21 N 

^^Ai„Ci<l'i(-Xi,z) + Y 5Z^*'"^*'^2(-a;i,2;) = 1 (4.3) 

m=l i=l m=l+l i=l 

which is completely in agreement with the previous normalization due to Riemann-Roch 
theorem. Indeed, the expression on the left-hand side is a meromorphic function on the 
spectral curve, it has g+2l—l poles and it has fixed values 1 at the points Pi, i.e. additional 
21 conditions. This function has to be equal to 1 identically. Now all the preliminaries are 
done and we turn to the construction. The Lax operator and its eigenvector normalized 



9 



by the condition ( [4 .31 ) are natural coordinates on the phase space and we take as 5L 
and 6C the corresponding cotangent vectors. We only introduce new variables fc"'^"', 
such that k"'^'^ = k — 2({z) + ({z — zi) + ({z — Z2) which is a well defined function on the 
spectral curve, new variables A"'^"' = \i — ({ri){{Bl^A\) + {Bl^Al)) and the corresponding 
expression for the Lax operator: 

Lr = + az - ^2)((5r^2) - 1) - - z,)m-^A\) + 1))+ 

(1 - A*2)$2(xi - (4.4) 

The index new is omitted in the following. We also introduce the differentials on the 
total space 6k and 6z. To relate it with the canonical differential on the spectral curve 
we notice that the total space is the fibration on the specific moduli space with fiber 
which is the symmetric power of the spectral curve. We also notice that the sub-varieties 
z = const are transversal to the fibers. So we could extend differentials on the fiber to the 
total space defining them zero along the sub- varieties z = const. It means that we choose 
the connection which is zero along such sub-varieties. For this choice of a connection we 
have Sz = dz, where we use the notation d for the extended differential (for more details 
see 



14, ra 



The main object of our concern is the two- form on the phase space 

1 ^ 

u; = - ^ respX< C* {SL + 5k) A 6C >)dz, (4.5) 

where C* is the eigenvector of L satisfying < C*C >= 1. 
Theorem 4 The two-form ui is 

25^5^(7,) A 5A;(7,). (4.6) 



Proof We first clarify the meaning of this formula, k is a well defined function on the 
curve and 2; is a multi-valued function. After the choice of a branch of the function z, its 
values at the points 7^ are functional on the space of spectral data, i.e. on the phase space 
of our system. The only thing to justify is that the corresponding cotangent vector SzIjs) 
does not depend on the choice of branch. Indeed, the difference between such functionals 
on different branches depends only on the modular parameter r of the base elliptic curve 
which is a priori fixed. The differential Q =< C*6L A 6C > dz is meromorphic on the 
spectral curve, the essential singularities reduce. The sum of the residues at the points 
Pi is the opposite of the sum of other residues, fl has poles at the points 7^. Using the 
condition < C*C >= 1 we obtain that the residues at these points are 

res-y^n = - < C*5LC > A5z{'js) = -5k{-fs) A Szi'js)- 

Other possible poles of Q are the branch points of the curve, at these points C* has poles 
due to the normalization, but the differential dz has simple zeroes at these points. ■ 

This proposition shows that our two-form is symplectic, i.e. it is closed and non- 
degenerate. It is non-degenerate because the number of Darboux coordinates is equal to 
the dimension of the phase space. 
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Theorem 5 The symplectic form lu admits the representation 

N 



cj = ^(5Ai A 5xi + SBf A SAl - SBf A 5A]^). (4.7) 

i=l 

This is proven by a straightforward calculation. We use the notation: 

Bi = {Bli, ■ ■ ■ , Bli, B^^, ■ • ■ , Bl). 
Let us decompose L,C as L = GLG^^, C = GC where G = Diag{e''^^^^^). Then 

a; = Kes,=oTr[6L A6h + C-\6L6C + [6h, L]a6C- 6hC A 6k)]. 

Here the expression G corresponds to the section of the direct image of the line bundle 
on the spectral curve to the base elliptic curve and Sh = G^^6G = diag{—5xiQ{z)). Using 

Tr{G-^[5h, L] A 5G) = -Tr{G-^6h A {6LC - C6k)), 

we obtain 

u = Kes,=oTr{2SL A Sh + C~^6L A 6G). 

The first summand gives Yli^^i ^ ^^i- Calculating the second summand we take into 
account the normalizing conditions ( [4.2|) 

N N 

Y.i^^o.Sc'; + 6A^^cf) = 0, J](4/S,„ + 54,5°) = 0, 
which demonstrates the proposition. ■ 

Theorem 6 The system of poles dynamics of the finite gap potentials of the Davey- 
Stewartson equation ^2.1\ ),...,( \2.13i ) is hamiltonian with respect to the symplectic form 
uj and the Gaudin Hamiltonian 

h = \Y. ^^(^r^i + B\^A^,) - avm''A\ + im^Ai - 1)+ 

i i 

^(i?r^^,)(i^r^i) "["(7"-+^^ (4.8) 

^ CT(r/)0r(Xj- - Xi) 

is equal to {Hi — H2)/2, defined in section 1 for the Gaudin elliptic system. 

Proof We substitute the flow defining the dynamics, i.e. the Lax equation, into the 
symplectic form. Our intention is to show that iQ^uj{X) = uj{X,dt) = dH{X). With the 
chosen coordinates the vector field is defined by the expressions: 

dt{Li-L,) = [Li,L2], dtk = 0; 
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dtC{t, P) = L^Cit, P) + /i(t, P)C(t, P) 
where P) is a scalar function. 

21 

2 



^g^cu = l^resp,(< C*{5L- 5k){L2 + ii{t,P))C > - < C*[L2,Li]6C >)dz 



i=l 

21 



^respMn{t,P)dz. (4.9) 



i=l 



The function /i has the same expansion at the marked points as the coefficient at t in the 
exponential part of the vector \E'. It means that at Pj,z = 1.../ the principal part of is 
equal to |/c and at the points Pj, z = Z + 1, 2/, ^ — 



1 ' 1 ^' 



4 ^ ^» 4 

i=l i=i+l 



; ^ 21 



^ resp^5 < C*L^C >dz + ^ ^^^p^^ < CU^C > dz 

i=l i=l+l 
1 1 

= -reSz=z2Tr5L'^dz — -reSz=ziTrSL'^dz. 



Substituting the expressions for the matrix elements of the Lax operator ( |4.4| ) into this 
formula we find 

reSz=z2TrL?dz — reSz=z^TrL?dz = 

i 

We have found the Hamiltonian up to a constant. To obtain the Hamiltonian Hi — H2 of 
the Gaudin system we have to subtract N((r]) and use the condition {B{^ A^) — (Bl^ A\) = 
2. We only need to note that the system (pr7|),...,( |2.12| ) is equivalent to the system of the 



hamiltonian dynamics with this Hamiltonian, and the equivalence is given by the formulas: 

al ^ exp(- [ (Ai + Xi)/2)ai; 



The "action-angle" variables can be constructed in the traditional way as in [|16| . We 
introduce the differentials dflj,j = 1, . . . , 2/ — 1 of the third kind with residues 1 and — 1 
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at the singular points Pj and P21 respectively. On the fundamental domain of our curve 
we can define the functions Ak{Q), k = g + 21 — 1 such that 



rQ rC 
AkiQ) = duji,i = 1,. . . ,g, Ag+i{Q) = 



dVLi.i = 1,...,2/-1. 



Now we introduce the quantities 



9+21-1 

0i = ^ Aii^s), Ik= <P kdz, Ig+i = resp^kdz. 

S = l ^"-i 
The symplectic form is 

UJ = 6(f)i A Sli. 

i=l 

Notice that the quantities Jj are integrals of motion. So, the obtained coordinates are the 
"action-angle" variables for our system. For the s/(2)-case the separated variables was 



also realized in 17 



5 Degenerations 



The construction of the Hitchin system on degenerate curves is well suited for the analysis 
of several degenerations of dynamical systems. For example the rational and trigonometric 
degeneration can be obtained by considering the elliptic curve with cuspidal and nodal 
singularities respectively instead of the regular elliptic curve. The normalization in both 
cases is the rational curve and the moduli of holomorphic bundles are parameterized only 
by the elements gluing the fibers over the singular points. For the Gaudin elliptic system 
such degenerations exist. The rational one is represented by the hamiltonian 



1 



(5.1) 



with the canonical symplectic form uj in the representation (O). This system involves 
a rational curve with two simple points and one double point with cuspidal singularity. 
It also could be solved by the algebraic-geometric direct and inverse problem as it was 
done in |12|, [|1^. The expressions for the \E' function could be found in [|13| (2.53). 
The "action-angle" variables for it were also constructed in JlSl. The author deals with 
another Hamiltonian but the separated variables for these systems coincide. There is 
another type of degenerations of our original system. It is the limit zi ^ Z2,'r] —>■ 0. This 
case is successfully taken care of by the direct and inverse problem as well. The FBA 



function can be found in |T3| (2.60). The solution provides the double periodic potentials 
for the Davey-Stewartson equation. The system can be specified by the Hamiltonian 



Hell =Ij2 ^^(^r^i + i^r^2) + x;(^r^D(5r^i)c(a;. 



(5.2) 
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and the standard symplectic form ([4.7|) . A priori this system is not integrable. Proceed- 
ing with the analogous construction for the "action-angle" variables, we obtain that their 
number is 2/(2iV— 2/-I-1) which is less than the dimension of the phase space. The same de- 
generacy can be observed by counting the number of Hitchin Hamiltonians. The Hamilto- 
nians of the original system can be represented as H^''- = res z=zg{{z — Zs)^^^Tr^'') , s = 1,2. 
When the two fixed points coincide, there is no more residues of the expressions Tr^'', 
and the only retrievable quantities are iff'' + H2'''. This case should be investigated with 
a more delicate geometrical analysis of the jet structure of the holomorphic bundle. 



Conclusion 

The main result of this paper is the explicit solution of the Gaudin elliptic system with 
spin and the uncovered connection between this system and the quasi periodic solutions 
of the Davey-Stewartson equation. For the methodological aspects we noticed the crucial 
correspondence between the Hitchin and finite gap descriptions of the dynamical systems. 
Our subsequent concern is to give a more general illustration of the equivalence of these 
two settings. The important technical query in this direction is the geometry of the mod- 
uli space of semi-stable bundles when we vary the base algebraic curve. It could clarify 
the Hitchin description of such a system as an elliptic degenerate Gaudin model. Such 



an attempt was undertaken and the case of the elliptic curve was analyzed in [|T9|. The 
construction crucially used the group structure on the elliptic curve and could not be 
translated directly to the general situation. In order to realize the geometrically estab- 
lished limit procedure a kind of canonical connection on the bundle over the moduli space 
of algebraic curves with moduli space of holomorphic bundles as the fiber is required. 
And the role of the Painleve connection in this context is of immediate concern to us. 
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